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The discrete static properties of a class of deformable double-well potential models are inves- 
tigated. The Peierls-stress potential of the models is explicitely given. Numerical analysis of the 
equation of motion reveal different soliton wavetrain profiles, most of which are periodic. Soliton 
wavetrains are also found analytically in terms of continuum nonlinear periodic wavefunction solu- 
tions then called periodons. The periodon stability in a lattice phonon bath is discussed. Looking 
at bifurcation phenomena and routes to chaos for a representative model of this class, the return 
■ map is derived in terms of a two-dimensional, two-parameter map. This map appears to be area 

^q' preserving, possesses three characteristic fixed points with one elliptic, and displays complex bifur- 

cation diagrams with hopf-like singularities. The routes to chaos also show complexities all due to 
the interplay of two control parameters. 
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I. INTRODUCTION 



The soliton concept (T]-|| is now well understood. However, in the past the subject casted great controversal, but 
also rich debates around two main themes: the soliton proper existence as real physical object and later, its stability 
namely against collisions with other objects of the same kinds(solitons, phonons) and with hosts(isotopic impurities, 
etc.). The concensus which emerges from these past activities involves the recognition of the universal role that soli- 
ton and in general, nonlinear waves play in various applications. In condensed matter, the soliton concept attracted 
attention for its intrinsic interest in the fundamental problem of phase transitions pj. In short, the main feature here 
is that two or more phases of materials may coexist separated by well-defined phase boundaries. As phase transition 
■ processes are dynamical by nature, the great theoretical challenge has been to find appropriate scenarios. The most 
popular one consists to manage at least one of the possible phases and next to follow its evolution among different 
• other phases. This scheme takes into account the nucleation, propagation and anihilation of phase boundaries. Dur- 
ing nucleation, the transition involves by turns an homogeneous, single-phase configuration and an inhomogeneous, 
multiphase configuration. The anihilation takes place in the opposite direction. Propagation involves motion of an 
already existing phase interface, which moves through the system increasing or decreasing volume of one phase at the 
, expense of the other. 

While this scenario applies in various theoretical contexts, the physical systems called in question are not always the 
same. For this reason, several one-dimensional(lD)models are used provided acceptable basic structural considera- 
tions. In the theory, each model bears on a specific one-site potential expected to provide the appropriate features for 
agiven physical system. The one-site potentials which have been widely studied are the </> 4 fffi, and sine-Gordon(sG) 
jApl. Both admit nonlinear solutions which are solitons. By its shape profile, the kink soliton has quite often been 
associated with phase boundary phenomena P,p|,[To[ . 

In practices, phase transitions are actually more complex and phase boundaries can largely depart from the sim- 
plest " single- kink" view-point pdjl. To this subject, relevant insights have been gained in some interesting results 
namely of Aubry fl2|| , Bak et al. |13|| and Pokrovsky jL4j, from which we learn that discreteness of the propagation 
medium is an intrinsic deal of the problem. Thus, lattice discreteness can by themselves determine the becomings of 
phase boundaries. This leads to relevant phenomena as "pinned interfaces" and "domain" configurations, revealed 
by experiments. In addition, the transient phases can display unexpected behaviours then indicating (phase and/or 
shape)memory-loss processes. These unpredictible behaviours are referred to as chaotic regimes, also pointed out and 
discussed at length [p~5| |l8[| . 

The recent advance on the soliton theory was directed toward new models with parametrized one-site potentials 
||l9| . Their advantageous features can be summarized in the fact that they reproduce the <fi 4 and sG while avoiding 
must of their shortcomings. This paper deals with a class of deformable double-well potcntials(DDWP) introduced 
very recently ||^,|2^,E3]. None of the versions of this class is equivalent to the familiar sinh-Gordon(shG) potential 
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j25[ , instead they must be regarded as distinct generalized versions of the shG model. These DDWP have several 
virtues including analytical tractability, in contrast with many others for which it is not always easy to get deeper 
in their explicit treatments. Moreover, their "double- well" shapes can be manipulated at one wish and hence allow 
fitting results to a desired context. To set these interesting features in a larger perspective, in the present work we 
will address to some questions related to lattice discreteness. 

The outlines of the work are the followings: In section II, we review their continuum dynamics and derive explicitely 
the Peierls stress potential (so-called Peierls-Nabarro, PN potential) by the Fourier-transform method |l^-|24|,|2(| . In 
section III We examine some configurations of the nonlinear solutions given by the numerical analyzis. An approxi- 
mate analytical solution is looked for in terms of periodons, i.e. periodic soliton wavetrains. The periodon stability 
in an intrinsic phonon bath is discussed. In section IV the return map of one representative model is derived in terms 
of a two-dimensional, two-parametcr(2D2P) map evolving in the phase space through distinct bifurcation scenarios 
in which alternate periodic and chaotic phases. Section V present concluding remarks. 



II. THE MODELS, KINKS AND PEIERLS-STRESS POTENTIAL 



Consider a one-dimensional(lD) Klein-Gordon(KG) system, that consists of a chain of harmonically coupled parti- 
cles with uniform mass M, each in the well bottom of a doubly-degenerate minimum potential. The dynamics of the 
model is described by an Hamiltonian: 



2 <Pn,t + 2 p 



'M+l 



(i) 



Where 4> n is the particle local displacement field, C a and u> the bare velocity and frequency respectively, and V(4>, fi) 
the DDWP that we assume of general form p3j: 
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The quantity a Q will be constant while fj, can vary, q and a are depending on the main deformability parameter /i 
and can take three distinct values i.e.: 





a = fx, q = Sinh~ 1 (/j) (3) 

For these values, V(6, ji) is a DDWP whose degenerate minima are continuously shifted but leaving unaffected 
the barrier height [f20|. 



ii ) 



a = Sinh q = Jt — = (4) 



For these second values, only the potential barrier will be affected while positions of the degenarate minima 
remain bounded at <p\ 2 = ±1 P2[| . 
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For these last values, the whole " double- well" shape of the potential changes by varying (i p3| . 

In the continuum limit, the generalized DDWP (||) can be treated analytically and yields the following single-kink 
solution: 



vt, u) = ±— tanh 
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The associate mass as well, in virtue of its particle-like properties, is analytically obtainable and reads: 

M fe ( M ) = ^L- [SinA(2a) - (1 + 2/i 2 )a] (8) 

This analytical expression of the generalized kink topological mass deserves instructive remarks concerning its defferent 
behaviours with repect to /x, for the three cases listed in i) — Hi): in the first, the kink mass displays a drastic decrease 
when increasing fj,, but remains always positive and non-zero. Similar behaviour prevails in the second case. In the 
third, by contrary, the kink mass behaves oppositely, increasing from negative to positive values through massless-kink 
shape for fi at about 0.76 Q. 

We now turn our attention to the response of the single kink obtained above to lattice discreteness effetcs. We can 
restrict ourselves to the static- kink regime, which is suitable for a "pinned- kink" configuration where kink pinnings are 
tough to result from forces provided by the discreteness of the propagation medium. These forces are often identified 
with the Peierls stresses governing structural changes in solid-state systems [[12-14 2flL T he generating potential is 



so-called Peierls-stress potential and has been the subject of a great recent interest[p7[, most often referred to as 
" Peierls-Nabarro" (PN) pinning potential. 

In the static regime, the total discrete system energy reduces to: 
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and the PN potential then follows by Fourier transforming G(n), i.e. calculating: 

u PN = Mci I™ dn iir) ex p( 2 ^ mn ) ( 10 ) 
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with: 
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It is easy to check that G{n) is an even function of the variable n. Therefore the Fourier harmonics follow from the 
class of integrals: 

dx c -2^n ff = ^ r ( 2 +; y ) r 3 ( 2 2 -^p-3g + . y[C0Sfe (x o )] (i2) 

[cosh(X) + co,sh(X )] 2 V2 T(2)sinh 3 / 2 (X ) s -1/2+*^ V o,\ y , 
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X = -= , cosh(X ) = 1 + 2/i 2 , Y = W2md(ji) (13) 

p in this relation is the Associate Legendre function here appropriately defined as: 

F 1 
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du[cosh(u) + cosh(X )]cosh(TrY)cos(Yu) 



While the complex- argument Gamma function T is such that: 



T(2 + iY)T(2-iY) = 



T(2) 



llfc=0 

with the product function evaluated by standard formula to give: 

II l +m 

k=0 

From now on, combining (^|)- (plf), we find: 
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Comparing the zero-harmonic U and the kink rest mass obtained previously in (g), it appears that this term repro- 
duces the kink proper energy. The next harmonics will then account for the lattice discreteness effects. The leading 
term U\ is dominant since the cosh argument is already almost insensitive to the second order in m. 



III. SOLITON WAVETRAINS AND PERIODONS 



When investigating the PN potential in section II, we assumed a continuum kink wave solution owing to the difficulty 
of proceeding with analytical treatment of the discrete motion equation. The best way of accessing discrete solutions 
is by numerical analysis. For this purpose, we will use the following nonlinear difference equation: 
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where d Q now sets the kink(or in general, soliton)length scale. We will only show the results for one representative 
among the DDWP above, that is |§. On fi gures 1 and figures 2 we displays some characteristic features emerging 
from the iterates in the phase space and of the wave shapes, respectively. It is manifest that the "wavetrain" feature 
prevails in the almost whole discrete regime. In particular, the discrete nonlinear regimes(Figs.2a, 2b) appear to 
be essentially dominated by " soliton-lattice" structures. These soliton-lattices solutions describe multi-kink states 
spreading along the chain axis with well defined periodicities. By varying the deformability parameter fj,, the "kink- 
lattice" structure collapses gradually and turns into almost harmonic periodic wavetrains via periodic structures that 
we can identify as weakly nonlinear, discrete modulational wavetrains(Fig.2c). 

Tough characteristic of the discrete nonlinear equation, the kink-lattice structure is not specific to the discrete phase. 
Indeed, solving the continuum limit of the nonlinear difference equation provided appropriate boundary conditions 
such as periodic boundary conditions, we may arrive at a continuum soliton wavetrain solution then called periodon. 
In principle, there are two different ways to check that the periodon describes lattice of kinks with well defined 
periodicity. The first is by direct integration of the motion equation ]28|-|3l|], leading to Elliptic functions which are 
indeed periodic. The second is by a naive and artificial method that consists in an algebraic sum of several(infinite) 
one-single kink wavefunctions p^ |. Both ways generate the same wavetrain profile and hence are equally valid. 
However, we choose to proceed via the first and find for the generalized version (pf): 
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Sn is the Jacobi Elliptic function of modulus A, K(X) the Jacobi Elliptic Integral of the first kind and K a is self- 
consistently given by the second-order polynomial: 



a 2 (l + ^ 2 ) 2 (l-A) 2 a 
64^ 4 + 8 



The period of the nonlinear solution (f20|) is: 
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To clearly see some interesting asymptotic behaviours of (|20|) , we consider the two limits A — > and A — > 1 . The first 
leads to: 



(s = x — t>i, /Lt) = ±— tanh 
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(24) 



This solution describes a breathing-like excitations, i.e. a kink with internal periodic quasi-linear oscillations. The 
effect of /j on the shape of this excitation mode is explicit from an expansion of (]24|) about small values of this 
paramter. The leading term yields: 
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For the second limit, expansions of the periodon wavefunction for small values of fi give back the single-kink solution 
(§• 

We would like to end this section with a short discussion on the stability of the periodon kink above against collisions 
in a phonon bath. Let us recast the periodon in the form: 



x — vt, u) = ± — tanh 1 [u(s)l, u(s) = B^Sni — =— | 



< A < 1 



(26) 



where the function u(s) is the part of the periodon wavefunction undergoing "shape dressing" upon collisions in the 
phonon bath. Thus we have to introduce the periodon-phonon ansatz as follows: 



u(s) = u{x) + t/j(x)exp(iu)t) 



The resulting eigenvalue problem is then: 
tp zz + [P(Q, A) — QXSn 2 (z)\il> — 0, z = 
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P(fi,A) = (l + fi 2 )(l + A) (28) 



This is just the Jacobian form of the (second-order) Lame equation 33 34]. The same equation is found for the <fi 4 
in the finite-support (or finit period) limit |33|. Its bound state spectrum possesses four non-zero frequency modes in 
addition to the translational mode. 



IV. BIFURCATIONS AND ROUTES TO CHAOS 



Modulated systems appear as the most frequent candidates when evoking systems that can be modelled within the 
present theory. These systems have long served as prototypes of materials with complex behaviours. The modulation 
periods are either commensurate or incommensurate with the basic lattice then providing them with their unique 
and rich physical properties. The transition between the two phases involve domain walls which are pinned soliton 
lattices. Properties of these soliton lattices are fundamental in determining the nature of the prevailing phase. It 
has been shown [[fl]-[l8| that the appearance of the soliton lattice can occur simultaneously with new phases which 
are essentially metastable or strongly unstable. Among them, lock-in and chaotic regimes correspond respectively 
to periodic domain-wall lattice and randomly pinned domain-wall lattice configurations. Theoretical attempts to 
understand these chaotic regimes have been carried for both <fi 4 and sG, in connection with physical phenonema as 
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Peierls instabilities. 

To a general viewpoint, the possibility of describing classical mechanical systems by discrete difference equations 
provides an excellent framework where to connect the well-understood periodic and chaotic behaviours of deterministic 
models with the equivalent behaviours in discrete nonlinear KG systems. The discrete logistic map is a sound 
ground for this purpose. According to existing approaches, the system will be said to exhibit chaos when its phase 
space trajectories show initially exponential sensitivity to small differences in initial conditions. This feature can be 
made quantitative by showing that the largest Lyapunov exponent is positive. We already know that for classical 
Hamiltonian dynamics with sympletic structures, the sum of all of the Lyapunov exponents is zero such that at least 
one is positive, unless all vanish as they due for integrable systems. By these means, chaotic instabilities are virtually 
generic for non integrable conservative Hamiltonian systems. In a past work |pq| , Bak and Jensen have considered the 
problem for the (j) 4 , showing that the theory could be formulated as an area-preserving 2D mapping. Their approaches 
cover a wide range of models involving those dealth with in the present context. Still, the competitions between the 
two main parameters(a and (i) in the DDWP are likely to generate new singularities in both bifurcation diagrams 
and chaotic windows and it is the aim of this section to point out some of the revelant ones. 
Starting, consider the 2D2P map: 
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M( QjAl ) can evolve among stable and unstable equilibra which are the fixed points of the map, and any set of points to 
form a trajectory on the map area will result in a periodic or chaotic orbit of one of these fixed points. As it sounds, 
whether the fixed points are stable or not is relevant in determining the local behabiours of the map area. Thus, a 
stable fixed point will be such that the orbits near it are moved even closer under the map, whereas nearby orbit 
points of the unstable fixed point will move aways by iteration. We can understand these equilibrium behaviours 
in the simple viewpoint of linear stability analysis, which consists to look at the transformation under M( a ^ as 
"gradient crossings" of neighbouring orbit points on the surface area consisting of the whole trajectory set of the map. 
Mathematically, we end with the matrix equation: 



(a.fj.) 



M n W n = W n 
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Where M n , then called Jordan matrix of the map, is the square Jacobian matrix of the gradients components at 
points W n = [u n ,v n ] on the map area. From now on, further iteration will result in conjugating the matrix M n 
such that for a given period-^ orbit of some fixed point W , the generator of the period-Z orbit points is M a . It 
turns out that properties of the transformation matrix will have direct incidence on the geometry of the map area. 
Indeed, since a matrix with determinant of absolute value different from one traduces change of norm, for the map 
this will correspond either to the contraction or to the expansion of the map area about the fixed point. Else, the 
map will be invariant under transformation and hence is area preserving. On the other hand, if all the eigenvalues 
of the matrix have magnitudes less(greater) than one, the associate fixed point will always attract (repel) the orbit 
points and therefore is and attractor(repeller). If one eigenvalue is less than one and the other is greater, the fixed 
point is a saddle: any perturbation of the orbit away from it may be magnified by further iteration leading to an 
unstable fixed point. The fixed point will be hyperbolic if none of the eigenvalues has magnitude one, and elliptic if 
the two eigenvalues lie on a unit circle one the map area. In theory, bifurcations of any of these fixed(or of any other 
periodic)points on the map can occur only for the parameter values for which an eigenvalue of the corresponding orbit 
set has absolute value one. Basing on this preliminary, we start by looking for the fixed points of the map M( a fi y 
They are three: 
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Sinh Sinh 
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It is noted that they coincide exactly with the fixed points of the ID invariant subspace of the map. 
matrix for the period-one orbit set is the 2D matrix: 



(31) 
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The norm of its determinant is one whatever the values of a and \i: therefore the map is always area preserving. Its 
eigenvalues at the fixed points W°, W l are respecively: 
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Where the signs =F stand for i=l, 2. While the fixed point (0, 0) is elliptic, for the two others the asymptotic behaviours 
are not so clear. We can circumvent this uncertainty by considering instead the traces of the matrix at those points. 
We obtain: 



Tu ,=2 
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In terms of these traces and according to our preliminary statments, the stability condition can now turn into the 



inequality: I T n 
hold for T+ 

we select the following parameter values: 



L i(a n) 1^ 2 ' Since M an< i a are au positive definite, it results that this condition on the trace will never 
hold for Tjt , and therefore W + is an always unstable fixed point. By contrary, | T~, ^ \ can be less than 2 provided 



a > 



y/l + fi 2 Sinh(2fi) 

In which case W~ is a stable fixed point. Otherwise, it is unstable. The marginal stability at 



yjl + n 2 Sinh(2n) ' 



(37) 



(38) 



marks a critical state at which the fixed point can bifurcates carrying the orbit either aways from it or closer by 
changing period. Since similar features hold for the fixed point (0, 0), everywhere below we will deal with trace. 
The stability of the next orbits, i.e. the period-two orbit and higher can from now be easily analyzed. Besides we are 
able to construct the recusion laws governing successive orbits around each fixed point. Indeed, writing the matrix 
equation for the transformation to a given 1-cycle orbit as W = M l W, we find by turns: 



rpQ c\l 
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2 i 
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(39) 



We learn from these recursion relations that the sequences of orbit cycles are indeed all periodic and therefore result 
in the map dynamics characterized by changes in the geometry of the KAM surfaces. As the phase-space trajectories 
in figures 1 illustrate, these per-iterate changes involve stretchings of the two axes as the separations between nearby 
points increase or decrease. The stretching factors in each of these two directions can be identified as the largest 
Lyapounov exponents of the orbit. For the fixed point (0,0), the recursion suggests that the norms of the largest 
Lyapunov exponents of successive orbits all converge to the same limit L — ln(2). Positive Lyapunov exponent 
traduces the fact that bifurcations of the fixed points will consist of period doublings(or halvings) as long as orbits 
remain asymptotically periodic, and can involve openings of chaotic windows along the bifurcation loop. By the same 
means we can estimate the convergence limits of the largest Lyapunov exponents of the fixed points . we find one 
single expression for the two non-zero fixed points @)- (|H|) i.e: 



L {l) = ln(2) + In 



1 



j (l + M 2 ) 
A* 2 



Sinh 2 {2n) 



(40) 



According to figure 3a and 3b , depending on the values of a and \x the parameter can change from positive to 
negative values and crosses a zero value without discontinuity. This behaviour reminds Hopf singularities. 
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To fix our minds as concerns the map sensitivity to the two parameters a(a — M M map) and — M a map) taken 
separately, we follow an arbitrary trajectory set from an input position for some range of values of each of these 
parameters. This leads to the bifurcation diagrams displayed on figure 4 and figures 5a-5c. Figures 4 corresponds 
to the a — M M map, with input u Q = 1CP 7 and for /j, = 0.5. Periods of the orbits seem to get halved converging 
to the critical value a — 1.685. Next, the first period doubling points out at a = 1.963, followed by a cascade of 
further period-doubling bifurcations with alternating chaotic and periodic attractor windows. Namely, the period-six 
attractor window is observed at about a = 2.165. Similar features also emerge for the \i — M a map. However, here 
emerges two main bifurcation scenarios shown on figures 5a and 5b, respectively. Figure 5a corresponds to values of 
a in the range a < l(in this particular case a = 0.99). The first crisis consists of period doubling at /i = 0.76 with 
the period-six chaotic attractor windows settling at fj, = 0.937. For a > 1, the prevailing bifurcation scenario is that 
shown on figure 5b(a = 1.99). On figure 5c we magnify the region of figure 5b that shows apparent crisis. Here also we 
note "period-doubling" cascades as well as chaotic phases marked by the random points. These complex behaviours 
suggest that the /x — M a map may show deeper sentitivity to parameter changes than the a — M^. 



V. CONCLUSION 



In summary, we have explored the essential discrete properties of a class of DDWP models. We have been able 
to derive analytical expressions of some important physical quantities as the PN potential, which makes the DDWP 
analytically tractable an therefore applicable to a large variety of discrete nonlinear systems. The multiplicity and 
richness of the distinct discrete wave amplitude solutions resulting from the numerical simulations, as well as the 
bifurcation diagrams showing the coexistence of period-doubling and period- halving cascades as well as chaotic phases, 
are of great importance for the possible large variety of domain-wall and phase regimes that they allow to explore 
within one single model. 
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VI. FIGURE CAPTIONS: 



Figure la: Phase-space trajectories for /j, <~ 1.5, a <~ 0.8. 
Figure lb: Phase-space trajectories for /i <~ 0.1, a <~ 0.8. 

Figure 2a: \i ^ 2 . Well-shaped discrete periodic soliton wavetrains given by numerical simulations. 

Figure 2b: fi ~ 0.5 . Discrete periodic wavetrain of narrow solitons. 

Figure 2c: /x ~ 0.01. Weakly nonlinear, discrete periodic envelope modes. 

Figure 3a: Largest Lyapunov exponent of the W fixed points as function of /i. 

Figure 3b: Largest Lyapunov exponent of the fixed points as function of a. 

Figure 4: Bifurcation diagram of the a — Mi a ^\ map. 

Figure 5a: Bifurcation diagram of the \x — M( a ^^ map for a < 1. 

Figure 5b: Bifurcation diagram of the /x — Mr a ^\ map for a > 1. 

Figure 5c: First critical region of the Bifurcation diagram of [i — M^ a ^y 
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